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Abstract 

We find that if the Higgs mass is close to its present experimental lower limit 
(100 GeV), Yukawa interactions in the quark-Higgs sector can make substantial 
contribution to the heavy quark MS masses. 
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1 Introduction 



Recent reports from ALEPH, DELPHI, and SLD |]J of measurements of the MS mass 
of the b-quark at the mass scale of the Z mass, combined with previous determinations 
at the scale of the b mass , show the running of this mass, and invite comparison 
with theory. So far, this has been done entirely in the context of QCD. In this paper 
we note that because the t-quark is so heavy its Yukawa coupling to the scalar sector 
(including the charged Goldstone bosons in Feynman-t'Hooft gauge) is comparable 
to the strong QCD coupling constant at the scale of Mz- Here we consider what role 
these nonQCD interactions may play in the analysis of these observations. 
II One-loop perturbative treatment 

We work in a truncated version of the standard model in which only the strong 
QCD gauge coupling (g 3 or g s ),the Yukawa couplings of the heavy t and b quarks, 
and the quartic self coupling of the scalar fields are kept nonzero. The weak gauge 
couplings (gi, #2) and all the other Yukawa couplings are taken to zero. (and nib "C m t 
is set to zero except when it is the object being determined). For physical values of 
these masses and coupling constants, these approximations are too crude for preci- 
sion calculations, but they provide dramatic simplification and hence are very useful 
for orientation for more heavy duty calculations, and in some cases, are sufficiently 
dominant to provide useful estimates. 

We compute the one-point and two-point functions of the model through one-loop 
order in MS renormalization and in on-shell momentum subtraction renormalization 
(OS MOM) in order to determine the counter terms required for these two different 
renormalization schemes. The relation between the perturbative pole mass and the 
MS mass follows from the pole condition 



We use m for generic renormalized mass,m for MS masses, and (sometimes) m* 
for perturbative pole mass. 
(1) t-quark MS mass 

Through two loop order for QCD and one loop for electroweak @@] , the 
relation for t-quark is 



= m* — m — S(m*) 
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+ ^f-2(-W) + \ ln 4 + ^r(l + 1- 4) - Tt4(3 + 3 In £)} (2) 
lb n z 2 4 16 7r^ 

Here N c = 3 and 

r = — , d(r) = --+ dx(2-x)ln(r 2 (l -x) + x 2 ) (3) 

m t 2 Jo 

For a perturbative treatment, we input to the right hand side the star quantities 
which are either measurable or expressible as combinations of measurable quanti- 
ties (We consider that the pole mass of a heavy quark can be determined up to an 
ambiguity of order Aqcd) 

m\ = 174 v* = 250 m\ = 4.9 (4) 

The Higgs mass is put in as one of a series of values ranging from 100 to 300 GeV. 
The coupling constants are fixed by the relations 

* = i = t,b, A = ^| (5) 

v 2v z 

In particular, y* = 0.9843 while y% is so small that we neglect its contribution to 
other quantities. These relations hold for all orders for bar and star quantities. 

The conventional definition of the MS mass is m t (m t ). With this definition, the 

2 

ln(^2-)'s in (g) are zero, and (g) reduces to 



CM A- C\ 

m t = m *{l - 5^(1) - p) 2 (15.37 - 1.04 n f ) 
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In Table 1. we give values of the nominal MS top quark ma,ss,m t (m t ) , computed 
with just the QCD contribution (a s ) and with the nonQCD contribution (y t , A) in- 
cluded. We see sensitivity to the Higgs mass. For the lightest Higgs considered, 
m h = 100, the nonQCD contribution is twice as large as and of opposite sign as the 
QCD contribution. For the heaviest Higgs, mh = 300, the nonQCD contribution is 
only a small correction to the QCD only result. 

(2) b-quark MS mass 

For the b-quark, the flavor independent QCD contribution to the relation is the 
same, but the electroweak contribution is different. 
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m b (fi) = m*{l - — (~ + In ^) 
7r 3 mi 
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_ ( _) 2((15 .37 - 1.04 B/ ) + (- - - + -n t ) ln£, + (- - - + 
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Now when we set // = m& to get the nominal MS b-quark mass, the QCD logs are 
again zero but in the electroweak contribution there remain large logarithms of 
These logs are not the result of running anything from fi = rrit to /x = w&. They arise 
from two features of the electroweak interactions which are not present in QCD. First 
is the existence of large mass splitting between members of a multiplet, in particular, 
t and b. The first In^- in (|7]) comes from the b-quark self-energy diagram in which a 
b-quark emits a virtual charged Goldstone boson and propagates as a t-quark until 
the boson is reabsorbed. In QCD gluons change color, not flavor, and quarks of same 
flavor, but different color have same mass. Second is the feature of Spontaneous 
Symmetry Breaking for mass generation and the concurrent existence of a scalar field 
which acquires a vev. When the scalar field is shifted to get into the proper vacuum, 
the vev propagates into the boson l-,2-,and 3-point functions and the fermion 2-point 
functions, and is the source of the second ln ^ and the ln ^ in (0). (For more detail, 
see the appendix). 

For the smaller Higgs masses considered (100,150 GeV), the perturbative mass 
shifts with n fixed at 4.9 are large enough that the condition m = m(fi = m) is badly 
violated, so we let jj, float and solve by iteration until the /x put into (^) is the same 
as the output rh(fi). The results are given in Table. 2. For Higgs mass of 300 GeV, 
the electroweak contribution is again a small correction to QCD contribution, but for 
lighter Higgs, including the electroweak contribution leads to substantially smaller 
MS b-quark mass than given by QCD contribution alone. For Higgs mass 100 GeV, 
the iteration is unstable and does not converge to any positive value. (The connection 
between a "light" Higgs and breakdown of MS perturbation theory has been noted 
before [Tl|). 

Ill Evolution of m&(/x) from \x equal to Mz 

From /i = 4.9 to 91.19 is a large change in scale, leading to large ln[— f-. We 

b 

therefor differentiate (|7|) to obtain the one-loop differential equation for the running 
of m&(/i) By keeping just the dominant one-loop QCD and electroweak Yukawa terms, 



4 



we can obtain a simple analytic solution. With t = ln( 



4.9 > 



Jt m b (t) - m b {j^(-*) + (— } = -m b A(t) 



m b (t) = m£exp(- / dt'A(t')) + const (9) 

•/to 

The initial condition 

rh(t ) = m 

fixes 

const = rh — m* 

Then we need the solutions of the one- loop RG equations for the QCD and electroweak 
Yukawa coupling constants. The one-loop QCD coupling constant is 



f 

l+9 2 so b l 

with 



'22 - -n, 



16tT 2 

and g 2 o = 2.599 fixed to make a s = .119 at fj, — Mz- 

For the electroweak Yukawa coupling constant, we start with the equation relating 
V and y* [J7| 



3/t =3/? (o - 2m -|, 



7T 3 U 2 

+i(-(^ c + f ) In 4 + f - d(r)) + ^(1) (11) 



Differentiating, the one-loop RG equation for y 2 is (y = y t = yt op to be distinguished 
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from t = In ^ 



d » 1 



^/ = ^((2^ + 3)(Z/ 2 - 12«^) (12) 

After substitution of g 2 (t) (|T0|)into ( pT2|) the resulting differential equation for y 2 {t) 
can be solved. 

V 2 (t) = y 2 {(l +1 t)U *jL(l + yt)* - (1 + TO)}" 1 (13) 

7 — 

with 



and the initial condition that y^(fx = m t ) is given by (O) fixes y\. With g^(t),yf(t) 
determined, the integral (|]) can be done numerically (t z = In ^ = 2.9237). 



Results are given in Table. 3. 

There is no entry in the table for = 100 GeV because there is no solution 
for that Mh in table. 2. The values of rhb{tz) obtained for Mh ranging from 150 to 
250 GeV fairly well span the range of experimental values |]] For Mh = 300 GeV the 
obtained value is probably high. Global fits to electroweak precision data |J disfavor 
Higgs mass greater than 200 GeV. For Higgs mass from 150 to 200 GeV, the values of 
rhb(Mz) in Table. 3 are consistent with the data while the values of fh^fhb) in Table. 2 
are not consistent with the generally accepted value (4.1 to 4.3 GeV). This leads us 
to question the provenance of the generally accepted value of ra&(ra&). It comes from 
QCD sum rule calculations, Lattice gauge calculations, and potential models. All 
done in context of pure QCD. Thus it is necessary to investigate the contributions 
of the electroweak Yukawa couplings to these calculations and see if they might then 
agree with the smaller values we have obtained for m&(wb) from the transformation 
(0), in the case of the smaller Higgs masses. 

The calculations just referred to are generally substantial undertakings without the 
added complication of the electroweak Yukawa contributions, and we do not propose 
to repeat and generalize any of them in this paper. However, we can identify one 
dominant term responsible for most of the electroweak contribution seen in Tables 
2 and 3, for the lighter Higgs masses. And to incorporate that term into a simple 
version of the QCD sum rule analysis is not a major undertaking. 

IV Electroweak contribution to QCD sum rule analysis. 

2 

Returning to (|7|) we see that for light Higgs the term j^"^f is almost the whole 
electroweak contribution. This term appears as an MS counter term for the quark 
selfenergy function. Its origin is in the elimination of the tadpoles introduced by 
the shift of the Higgs field to do perturbation theory about the stable vacuum. See 
Appendix A and 0. In the QCD sum rule calculations, one computes moments of 
the dispersion representation of the vacuum polarization. 




and 



m t {tz) = mle F + m b - ml 



(14) 




dxe*** < 0|T(i M (x)^(0))|0 >= {q^ - gV)(CV(g 2 ) + . . . 



(15) 



6 



cm = -91 r d3 J!^ (16) 

7T ./4m 2 S(S + Q 2 ) 

If 3?C(s) 1 /• , R{s) 

M » = «J dS ^=V2^J ^ (17) 

The QCD Feynman diagrams we include are the basic one-and two- loop vacuum 
polarization diagrams (contributions to the moments available in the literature ||) 
and two diagrams with a single — mass insertion in one or the other of the quark 
propagator lines in the first diagram . The new contribution to QC(s) is 



2N r 5 m m 4 



7T m S 2 Jl 4m^ 



with 



This gives calculated 



QSC{s) 



m \bn z r A mi 



M n = + + 8M^ ] (20) 



to be matched against the integral of R(s) evaluated from experimental data on 
the upsilon resonance plus continuum. 

In the purely QCD analysis, the convergence of the perturbation series is accept- 
able up to fairly high n moments. 

Recently Jamin and Pich || have done an exhaustive treatment of the purely 
QCD analysis. They found that the optimal choice of moments to employ was n 
from 7 to 15. (Smaller n is sensitive to poorly known continuum threshold, and for 
larger n perturbation theory gets worse and nonperturbative contributions come into 
play) Furthermore the value of the mass obtained is rather stable in this range. But 
precisely because the dominant MS electroweak term is so large,the perturbative 
treatment fails rapidly as the order n of the moment increases. We give the results 
for the moment A4\ in table. 4. For the integral over R in (|i~7|)we used PDG data for 



six upsilon resonances and asymptotic form Of R for continuum starting at 11.1 GeV. 
There is no entry in Table. 4 for Higgs mass 100 GeV. For this value of Higgs mass, the 
perturbative calculation gives an unphysical negative value for M. x So again we see 
the breakdown of the MS electroweak perturbation series for this "light" Higgs. For 
the other values of M/, from 150 to 300 GeV, the values of fh^frib) obtained from this 
crude implementation of the QCD sum rule, augmented by electroweak contribution, 
look at least recognizably like the results in Table. 2, obtained from the transformation 
equation (0). 
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V. Conclusions 

Although none of the numerical results presented here are definitive - the approx- 
imations have been too many and too crude, a consistent pattern has emerged, with 
nonnegligible electroweak contributions and values of m^ifh^) substantially smaller 
than currently accepted, and furthermore, sensitive to the mass of the Higgs boson. 
We have started from the position that we have a pretty good knowledge of the per- 
turbative pole mass of a heavy quark from threshold behavior,potential models, and 
lattice calculations. We have chosen 4.9 GeV (and don't protest if someone prefers 
4.8 or 5.0). We have then used the perturbative definition of the pole mass and MS 
perturbation theory to convert to the nominal m& mass rhb{rhb)- We pointed out some 
features present in the electroweak sector which do not appear in QCD. And we found 
values of m& which depend on the Higgs mass, and for range of Higgs masses favored 
by global fits to precision electroweak data, are significantly smaller than currently 
accepted value. At this point in the argument we are in conflict with one or the other 
of two accepted positions: Mh should be less than 200 GeV and mj is 4.1 to 4.3. We 
pointed out that this value for rhb is not directly a measured number, but depends 
on the analysis used to fit experimental data. In particular, the analyses which led 
to this value are all purely QCD. We then made an (admittedly crude) analysis of 
what value might come out of the QCD sum rule approach if one included the elec- 
troweak contribution. The results were gratifyingly similar to those obtained from 
the conversion from pole mass. We also carried out the evolution Of rhb from fi = rhb 
to /i = Mz, and found results consistent with the recently reported experimental 
results. There is a possible problem here. Just as in the case of rhb(fhb), the analysis 
presented in the experimental papers was purely MS QCD. I have not found a quick 
way to incorporate the electroweak contribution in this analysis, but I refer to a paper 
by Rodrigo |L0| which is aware of the possibility of important electroweak Yukawa 
contributions, but claims that they can be suppressed by taking appropriate ratios 
of ratios of the experimental data. Then our analysis of the evolution including elec- 
troweak contribution would remain consistent with the experimental results. Finally, 
the dominant large electroweak contribution we have exploited, 5( v , is specific to the 
MS scheme. As explained in the Appendix, in a theory with a global symmetry and 
Goldstone theorem,^* may simply be taken to be zero in a MOM scheme. There 
will still be a strong y t Yukawa coupling connecting the two members of the t,b- 
doublet, which suggests that even with 5Q = 0,the electroweak Yukawa interactions 
may make a contribution to the QCD sum rule analyses. From the point of view of 
a MOM scheme, the various m^'s are a nonissue, and the QCD sum rules become 
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a more general testing ground for the whole theory, QCD plus electroweak . In the 
absence of one dominant contribution, it will be necessary to do a complete pertur- 
bative calculation of the electroweak contribution to the QCD moment functions, to 
see how it works out. 



A details of renormalization schemes 

Since the 2N c y 2 (^) term is the primary source of our results which differ from the 

h 

received wisdom,we explain here in some detail how it arises and is propagated in the 



MS two-point and three-point functions 

The original Lagrangian for the Higgs-Goldstone sector is 

C = (d<S> B )(d<S> B ) - /4$L$b - Xb(^b) 2 (21) 

with 

For //q < we are in the broken symmetry phase and < Hb > is not zero in 
the stable vacuum. In order to be perturbing about the true ground state, rewrite 
everything in terms of the shifted field which has zero vev in the true ground state. 

<h B >=V B , h B = V B + h B , <h B >=0 (23) 

Renormalize 

V B = Jzlv (24) 

^B = ^m 2 ^ 2 ' ^ B = 

The quartic scalar vertex renormalization is Z 4 = Z\Z\ 

The complete all orders vev V is determined by the requirement that the vev of 
h be zero, order by order in the loop expansion. Our notation for the order by order 
expansion of V is 

& = l + f =1 + (C,-1) =1 + 5C K } 

We use the notation ( v rather than Z v , to emphasize that (]25|) is not the renormal- 
ization of any bare field or bare parameter in the Lagrangian (pip. It is the loop 
expansion of the vev of the renormalized field (and is UV finite in the truncated 



9 



theory considered here, but both UV divergent and gauge variant when the gauge 
couplings are turned on). 

The field shift ( |2"3"D generates in the Lagrangian ( |2"1~D terms linear in the shifted 
field. 

^linear = — V (Z '^Z '^2 /i 2 + Z±XV 2 )h 

= -v(fi 2 + Xv 2 )h - [((£, - 1)(^ 2 + Xv 2 )) - v{{Z+Zp - l)/i 2 + (Z 4 ( 2 - \)Xv 2 )\ h 

(26) 

Treated as a perturbing interaction, the first term in (|26| ) will generate tree tadpole 
diagrams. We eliminate them by fixing the tree level vev of the unshifted field as 

v 2 = Qx 2 < 0) (27) 

Then inspection of the terms in the Lagrangian which are quadratic in the fields after 
the shift, leads to identification of the masses, 

Ml = 0, M 2 h = M 2 = 2Xv 2 (28) 

Although this is a relation between quantities identified at tree level,in our treatment 
there are no higher order corrections. In an OS MOM scheme this mass is the physical 
Higgs mass (modulo problems with unstable particles) (fTTH). 

The field shift also generates in (21) terms trilinear in the scalar fields. There 
are also trilinear Yukawa couplings between the quarks and the scalars, which have 
been the subject of the bulk of this paper. These terms generate one-loop tadpole 
diagrams (fig. 1). After the condition (p7|) has been applied, ( p6|) is still a term of 
the Lagrangian linear in the shifted field, but the coefficients are formally of one-loop 
order, so the tadpole trees generated by it now act as counter terms to the one-loop 
tadpole diagrams (fig.l). 

4Lr = v{Z+Z? - Z 4 ( 2 )^h = v5f, 2 h (29) 
The elimination of the one-loop tadpoles fixes this particular combination of 5Z's and 

bix 2 = 3X{A + A M ) - 2N c y 2 A m , (m = m t ) (30) 

A is the standard tadpole loop integral, dimensionally regularized. 

^-i^^^^^f- 11 (31) 

2 

A e = 7£ + In 47r. e = 4 - d 

e 
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In MS all of the SZ's are only of form constant times A e and cannot cancel log and 
constant pieces of combined (|30|) , (|3~TD . Then the condition that the vev of the shifted 
field vanish in one-loop order fixes 



16tt 2 



, M 2 , 9 m 2 , , m? „ 

S^l-to )- 2*^(1 -In ^) 



(32) 



The field shift introduces V, and hence 5( v , into the quadratic terms in ( 21f) as 



well as generating the linear terms (pq). This is the origin of the large electroweak 
contribution to the transformation from OS to MS for Mh < m t . 

In OS MOM renormalization, any of the Z's, as well as 8( v . can contribute to finite 
part of (|30|). In addition to the condition of the vanishing of the vev of the shifted 
Higgs field, we have to specify a number of renormalization conditions equal to the 
number of bare parameters Xb) an d field multiplets (<&b) in the Lagrangian. We 
could take those to be the mass shell conditions (pole and residue) 

5^(0) = 0, SUM 2 ) = 0, £ h (M 2 ) = 0. (33) 



(For MS these conditions have zero replaced by finite, i.e. MS counter terms remove 
all A £ 's) In our notations, the two-point counter terms appear in combinations 



D h \q 2 ) = q 2 -M 2 - S^(g 2 ) + (Z, - l)g 2 + 5^ - (Z 4 ( 2 - 1)M 



The same 5fi 2 (|30| ) satisfies both < h >= and £</>(0) = 0. This is a manifestation 
of Goldstone's theorem and is a check on our renormalization procedure. From these 
conditions one can determine Z^, Z^and the product Z±(%. For MS , 5( v is separately 
determined by the vanishing vev condition, so all the Z's are determined (and it is 
verified that 5£ v is finite). For OS MOM scheme another renormalization condition 
is required to separate Z4 from ( v . For this, we can choose h(f)+<fr- vertex. This is 
related to the above two-point functions by a Ward Identity. 

VT hM _(0,q) = D^(q 2 ) - D h \q 2 ) (35) 

Separate this into tree and one-loop and also write separately integrals from Feynman 
diagrams and counter terms. 

(v + 5V)(2Xv + f FD + 2Xv(Z A ( v - 1)) = 

q 2 _ E FD + {Z<j) _ 1)q 2 + Sfx 2 _ q 2 + M 2 + J^FD _ _ 1)g 2 _ ^2 + (^2 _ ^ 

(36) 
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v2Xv = M 2 

vT FD (0, q ) = -Z FD ( q 2 ) + E FD (q 2 ) (37) 
SV2Xv + v2Xv(SZ 4 + 5( v ) = + (5Z 4 + 25( v )2\v 2 

The last line is an identity, The Ward Identity is satisfied for any values of SZ 4 
and 5( v . We complete the the renormalization conditions for the OS MOM scheme 
by the condition 

vT h4>+(t> _(0,q)\ q 2 =M2 = D^ l (M 2 ) (38) 

This condition requires SQ to be zero (V* = v*) and Z 4 is now fixed by (^) with the 
Ward Identity then implying ([38]). 
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Figure 1 



Table 1: t-quark MS mass mt(mt) 
M h QCD(l-L) QCD(2-L) QCD+EW 



100 


166 


164 


185 


150 


166 


164 


174 


200 


166 


164 


170 


250 


166 


164 


168.6 


300 


166 


164 


168 



Table 2: b-quark MS mass mb(mb) 
M h QCD(l-L) QCD(2-L) QCD+EW 



100 


4.46 


4.19 




150 


4.46 


4.19 


2.52 


200 


4.46 


4.19 


3.60 


250 


4.46 


4.19 


4.11 


300 


4.46 


4.19 


4.47 
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Table 3: b-quark MS mass at Mz, m\,{Mz) 
M h m t F m b (M z ) 
100 - 



150 


2.52 


-.665 


.0528 


2.27 


200 


3.60 


-.308 


.1679 


2.84 


250 


4.11 


-.196 


.2133 


3.17 


300 


4.47 


-.0918 


.2343 


3.45 



Table 4: b-quark MS mass mj, from sum rule 

M h QCD QCD+EW 

100 4.46 

150 4,46 2.86 

200 4.46 3.63 

250 4.46 3.98 

300 4.46 4.19 
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